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Abstract 

The spatial distribution of unvisited/persistent sites in d = 1 A+A — > model 

is studied numerically. Over length scales smaller than a cut-off £(t) ~ t z , 

the set of unvisited sites is found to be a fractal. The fractal dimension 

dj, dynamical exponent z and persistence exponent 9 are related through 

z(l — df) = 9. The observed values of df and z are found to be sensitive to 

the initial density of particles. We argue that this may be due to the existence 

of two competing length scales, and discuss the possibility of a crossover at 

late times. 
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Persistence properties of spatially extended systems undergoing time evolution has at- 
tracted a lot of attention of late. Generically there is a stochastic field </>(x, t) at each lattice 
site x, which evolves with time t through interactions with other (usually nearest neighbour) 
sites. One quantity of interest in the present context is the persistence probability at time 
t, which is defined as the fraction P(t) of sites in which the stochastic field </>(x, t) did not 
change sign in the time interval [0,t\. In a large number of cases, it is found that P(t) ~ t~ e 
|J. The new exponent 9, called persistence exponent, is, in general, unrelated to other 
known static and dynamic exponents. 

The non-trivial nature of 9 can be attributed to the interactions between neighbouring 
sites, which makes the effective stochastic process at any single site non-Markovian. Suppose 
0(x, t) flips sign at time t. This event will increase the chance of neighbouring sites also 
flipping sign at subsequent times t' > t. This leads to the growth of spatial correlations in the 
system, which die out with increasing separation, on account of the statistical independence 
of distant flips. The non-Markovian nature of the process is more directly captured in 
these spatial correlations. Although much effort has been expended for calculation of the 
persistence exponent 9 by exact |2| and approximate || methods, little has been done to 
investigate the associated spatial correlations in the process. In this paper, we undertake 
such a study in d — 1, where the correlations are expected to be most pronounced. 

In one dimension, the zeroes of the stochastic field can be viewed as a set of particles, 
moving about in the lattice, annihilating each other when two of them meet. When a particle 
moves across a lattice site for the first time, the field there flips sign, and the site becomes 
non-persistent. If each particle is assumed to perform purely diffusive motion, this reduces to 
the well-known reaction-diffusion model A + A — > Q , with appropriate initial conditions. 
The simplest case is random initial distribution of particles, with average density n , for 
which P(t) ~ t~ e with 5 = 3/8 @, independent of n ||. We investigate spatial correlations 
in persistence for this simple model. We start with the two-point correlator C(r,t), which 
is defined as the probability that site x+r is persistent, given that site x is persistent 
(averaged over x). We define p(x, t) as the density of persistent sites: ie., p(x, t) = 1 if site 
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x is persistent at time t, and otherwise. Then, with the previous definition, the expression 
for C (r, t) is as follows. 

C(r, t) =< p(x, t) p(x, t)p(x + r, t) > (1) 

where the brackets denote average over the entire lattice and < p(x, t) >= P(t). 

Our main results are the following. Strong spatial correlations exist in the distribution 
of persistent sites, with a cut-off length scale separating correlated and uncorrelated 
regions. At late times t (ie., in the scaling regime), this length scale grows as a power of 
time: ~ t z , where z is the dynamical exponent in this context ||. In the correlated 
region r <C £(£), the correlator shows a power law decay with distance: C(r,t) ~ r~ a . The 
scale-invariant behaviour, indicative of strong correlations, shows that the set of persistent 
sites is a self-similar fractal with dimension df = 1 — a. By consistency, the exponents 
are related as za = 9. We have analyzed the fractal structure by box-counting method 
also. Careful measurements of the exponents over several decades of MC time show that the 
observed values of a and z change with initial density n while satisfying the above scaling 
relation. 

We did our numerical simulation on a 1-d lattice of size N = 10 5 , with periodic boundary 
conditions. Particles are initially distributed at random on the lattice with average density 
n , and their positions are sequentially updated — each particle was made to move one step 
in either direction with equal probability (D = 1/2). Whenever such a move resulted in two 
particles occupying the same position, both are removed from the lattice before moving to 
the next particle. The starting density of persistent sites is P(0) = 1 — n , and a persistent 
site becomes non-persistent when it is occupied by a particle for the first time. The time 
evolution is done up to 10 5 Monte-Carlo steps (1 MC step is counted after all the particles in 
the lattice were touched once). These time and lattice scales are the largest possible within 
our computational resources. We repeated our simulations for a few values of starting density 
n . The results were averaged over 50 different initial realisations. 

For distances r ^> 1 and late times t, we find that C(r,t) ~ r~ a for r <C In 
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the other extreme of large separations, ie., r 3> £(t), the sites are uncorrelated so that 
C(r,t) = P(t) ~ t~ 6 , independent of r. Thus £(t) is the correlation length for persistence, 
and consistency demands C,(t)~ a ~ This implies a power- law divergence: ~ t z with 
a dynamical exponent z related to a and 9 through the scaling relation 



The observed behaviour of C (r, t) can be summarised in the following dynamic scaling 
form. 



with the scaling function f(x) ~ x~ a as x <C 1 and f(x) ~ 1 for x ^> 1. In Fig.l, the 
scaling function /(x) = C(r,t)/P(t) is plotted against the scaled distance x = r/t z for two 
values of time separated by a decade. The initial density is no = 0.5. Excellent data collapse 
is obtained for z = 1/2, and the measured value of the spatial exponent a ~ 3/4 is entirely 
in accordance with the scaling relation. 

The observed power-law decay of C(r,t) with r has a wider significance, apart from 
showing the strong spatial correlations in the distribution. It implies that, over length 
scales not too large, the underlying structure is a self-similar fractal. This is most easily 
seen with the 'box-counting' procedure f7j. We divide the enire lattice into boxes of size I, at 
time t. After discarding 'empty' boxes, ie., those which contain not even a single persistent 
site, let M(l, t) be the average number of persistent sites in a box of length I. This quantity 
is related to C(r,t) through M(l,t) = f C(r,t)dr. Substituting the scaling form Eq. |3] for 
C(r, t), one finds 



za = 6 



(2) 



C(r,t) = P{t)f(r/£{t)) 



(3) 



M(l,t) ~ / 



I « Z(t) 



(4) 



M(l,t)=lP(t) l^>£(t) 



(5) 



which can be summarised in the scaling form 
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M(l,t) = lP(t)h(l/Z(t)) (6) 

with the scaling function h(x) ~ x~ a for x <C 1 and ~ 1 for x S> 1. We see that 
over small enough length-scales / <C the set of persistent sites form a self-similar fractal 
with fractal dimension df = 1 — a, with a crossover to homogeneous behaviour at larger 
length scales. This crossover is illustrated in Fig. 2, where we have M(l,t) (measured from 
box-counting) plotted against the box size I for three values of time. The initial density here 
is no = 0.5, and we find df ~ 0.25 in agreement with our result from study of the two-point 
correlation C(r, t). 

In Fig. 3, we compare the results from box-counting for different starting densities. For 
n = 0.2, we see that the fractal region appears much later compared to higher values. 
This is presumably due to the large inter-particle separation at t — 0, and the consequent 
delay in reaching the scaling regime. For higher densities, the fractal dimension is seen to 
decrease continously with n , approaching zero in the limit n — > 1. We notice that although 
£(t) ~ 10 3 in terms of the lattice spacing, it is still much less than the lattice size N, so as 
to rule out finite-size effects. 

In Fig. 4, we plot the scaling function h{rj) = M(l,t)/lP(t) against the scaling variable 
rj = l/t z for two values of time separated by a decade. We have displayed results for n = 0.8 
and 0.95. For n = 0.8 the best data collapse is obtained with z ~ 0.45, wheras for n = 0.95, 
the corresponding value is z ~ 0.39. The exponent a, measured from the small argument 
divergence of h(r)), also shows similar changes. 

In Table I, we have summarised our exponent values for four initial densities. All mea- 
surements were made using the data for the mass-distribution M(l,t) rather than the cor- 
relator C (r, t) on account of lesser statistical fluctuations. For the dynamical exponent z, 
we chose the value which gave the best collapse of data under dynamic scaling. Although it 
is difficult to measure the exponent very accurately using this method, we have verified by 
visual inspection that the error involved is less than the reported changes in the exponent 
values at least by a factor of two. We have omitted the case n = 0.2 because no single value 
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of z was found to give good scaling behaviour in the time range studied. 

A more direct way to measure the dynamical exponent z is to compute the average 
separation L(t) between persistent sites. If the spatial distribution were uniform over all 
length scales, this quantity would be simply L(t) ~ Pit)' 1 . Since this is not the case, we 
have to proceed more carefully. We define n(k, t) to be the number of nearest neigbour 
pairs of persistent sites at time t with separation k. By definition, f k n(k,t) = NP(t) and 
J k kn(k,t) = N. The average separation L(t) = iV _1 J k k 2 n(k,t) and we expect L(t) ~ t z . 
We computed L(t) numerically by simulating 100 lattices of size N = 10 5 upto 10 5 time steps, 
for each no- In Fig. 5, we display the results for the running exponent z e ff = d(logL)/d(logb). 
The results are seen to be fully supportive of our earlier conclusions. 

Our numerical results are strongly suggestive of non-universal behaviour of exponents 
a and z. The non-universal exponent values have been observed to be valid over at least 
three decades of MC time (upto 10 5 time steps) . We note that there are two length scales at 
work here. For low n , the dynamics is dominated by diffusive motion of isolated particles, 
'eating into' clusters of persistent sites. Due to annihilation, their average density decays as 
n(t) = {SixDt)' 1 ^ 2 || and hence the average separation is the diffusive scale £o(t) ~ t x l 2 . 
On the other hand, for hq — > 1, the initial separation of persistent sites ~ 1/(1 — no) 3> 1. 
The short time behaviour is now dominated by persistent — > non-persistent conversion of 
isolated sites, with characteristic length scale C p (t) ~ t 3//§ . It is possible that the observed 
non-universal behaviour results from competition between these two scales. According to 
this picture, one should see a crossover to diffusion dominated regime at later times, but 
we are yet to see any signature of that. Further numerical work, at least a few orders of 
magnitude greater than what is reported here, would be required to establish conclusively 
the possibility of a temporal crossover. 

To conclude, we have discovered strong, power-law correlations in the spatial distribution 
of persistent sites in one-dimensional A + A — > model. The correlation length £ (t) exhibits 
an algebraic divergence with time. In the correlated region, the set of persistent sites form 
a self-similar fractal, while over larger length scales, the distribution is homogeneous. These 
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features are not specific to this model or dimension. We have observed identical features 
in kinetic Ising model in d — 1 and 2 ||, showing that this is a general phenomenon in 
the context of persistence. The interesting aspect of the present model is that the fractal 
dimension was found to be sensitive to the starting density of particles. 

We thank M. Muthukumar for discussions and G. I. Menon for a critical reading of the 
manuscript and suggestions. 
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TABLES 



no a z 

0.50 0.7342(8) 0.50 

0.80 0.8294(5) 0.45 

0.95 0.9517(3) 0.39 



TABLE I. Observed values of exponent a as measured from box-counting method (details in 
text), for four values of initial density hq. The quoted value of dynamical exponent z is the one 
which gave the best data collapse over three decades of time, t = 10 3 ,10 4 and 10 5 . The fractal 
dimension df = 1 — a. 
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FIGURE CAPTIONS 



FIG. I: The scaling function for two-point correlation f(x) = C(r, t)/P(t) plotted against 
the scaling variable x = r/t z on log-scale for two values of time t = 10 3 and 10 4 . The starting 
density of particles is n = 0.5. The data for different times are seen to collapse into the 
same curve if scaling is done with z = 0.50. The observed a ~ 0.75 is in agreement with the 
proposed scaling relation (2). 

FIG. II: The average number of persistent sites M(l,t) in a box of size / at time t is 
plotted against the box size I for t — 10 3 , 10 4 and 10 5 . The initial density of particles is 
n = 0.5. The crossover from fractal (dimension df ~ 1/4) to homogeneous (df — d — 1) 
distribution is clear from the figure. 

FIG. Ill: Same as Fig. 2, for four starting densities n = 0.2, 0.5, 0.8 and 0.95. All plots 
correspond to t — 10 4 . For n = 0.2, the fractal region is reached late, but the asymptotic 
value is seen to be the same as that for n = 0.5. For higher n , df decreases continously, 
approaching zero in the limit n — > 1. 

FIG. IV: The scaling function for the mass-distribution h(rf) = M(l,t)/lP(t) plotted 
against the scaling variable rj = r/t z for two values of time t = 10 4 , 10 5 and two starting 
densities no = 0.8 and 0.95. The observed data collapse has been obtained with z = 
0.45(n = 0.8) and z = 0.39(n = 0.95). The corresponding values for a are ~ 0.83 and 
0.95. For comparison, a straight line with slope 0.75 is also shown. 

FIG. V: The running exponents z e /f for four starting densities is plotted against 1/logi. 
These results have been averaged over 100 starting configurations. 
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